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We study the effect in geometrically frustrated antiferromagnets of weak, random variations in the 
strength of exchange interactions. Without disorder the simplest classical models for these systems 
have macroscopically degenerate ground states, and this degeneracy may prevent ordering at any 
temperature. Weak exchange randomness favours a small subset of these ground states and induces 
a spin-glass transition at an ordering temperature determined by the amplitude of modulations in 
interaction strength. We use the replica approach to formulate a theory for this transition, showing 
that it falls into the same universality class as conventional spin-glass transitions. In addition, 
we show that a model with a low concentration of defect bonds can be mapped onto a system of 
randomly located pseudospins that have dipolar effective interactions. We also present detailed 
results from Monte Carlo simulations of the classical Heisenberg antiferromagnet on the pyrochlore 
lattice with weak randomness in nearest neighbour exchange. 
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I. INTRODUCTION 

Frustration refers to competition between few-body in- 
teractions which hinders simple macroscopic long-range 
ordering. In many systems it involves competition be- 
tween ferromagnetic and antiferromagnetic interactions. 
However, frustration can also occur in systems with 
purely antiferromagnetic interactions. Geometrically 
frustrated antiferromagnets^ constitute a large class of 
materials in which the frustration has a purely structural 
origin and gives rise to highly degenerate ground states. 
In some instances a consequence of this degeneracy is 
that the system has no ordered, low-temperature phase, 
instead remaining in the paramagnetic phase down to 
zero temperature. 

It is well known that disorder in the form of 
quenched random-signed few-body interactions leads in 
high enough dimensions to a spin glass phase separated 
from the high temperature paramagnet by a true, if 
unusual, phase transition!^ It has been a longstanding 
question whether the addition of random interactions on 
top of those of a geometrically frustrated antiferromag- 
net could also lead to a true spin glass phase at low 
temperatures;^ and if so, whether and to what extent 
the transition and low temperature phase are similar to 
those of conventional spin glasses. 

In fact, spin glass like freezing has been observed in 
a number of geometrically frustrated magnets. These 
materials characteristically have a Curie- Weiss constant 
(^cw of magnitude much greater than the freezing tem- 
perature Tf. Some examples are SrCrgGa40i9 (6cw — 
~500K,Tf ~ AK)!^-! Y2M02O7 {Ocw - -200A', Tf ~ 
22K)^-22. and Zni_^Cd:rCr2 04 for x ~ 0.05 {Ocw ^ 
-390, Tf ~ 12K),i^ Typical observed features include 
differences between field-cooled and zero-field-cooled sus- 
ceptibilities below Tf, and in some cases an increase in 
non-linear susceptibility Xni close to Tf. In particular, 
the existence of a sharp spin glass transition has been 
rather clearly established through detailed experiments 



on Y2Mo207,a>a 

The reason for the observed freezing has long been a 
puzzle. On one hand, it has been established that sim- 
ple models without disorder do not show freezingji^ On 
the other hand, samples exhibiting spin glass order ei- 
ther contain little structural disorder— that could be in- 
voked to explain the transition, or have a transition tem- 
perature that does not correlate straightforwardly with 
the level of the identified form of disorder. Indeed, in 
SrCr8-xGa4-|_xOi9 the transition temperature decreases 
with increasing disorder, as represented by the composi- 
tion 

A possible origin for a low temperature spin glass phase 
in frustrated magnets is suggested by recent experiments 
that show the importance of random strains in the sam- 
ples. Such strains, via magneto-elastic coupling^ gen- 
erate randomness in the strength of antiferromagnetic 
exchange and hence may account for a spin glass phase 
at low temperatures. In the material Y2M02O7, disorder 
in Mo — Mo distances has been detected using XAFSi^S 
Separately, in Zni_2;CdxCr204 disorder can be intro- 
duced in a controlled fashion by varying the composition 
X. Since Zn^+ and Cd^"*" have different ionic radii, this 
non-magnetic disorder is expected to introduce random 
strains. Moreover, the fact that the undoped material 
{x — 0) has a low-temperature phase transition at which 
a frustration-relieving lattice distortion and Ncel order 
appear together— suggests there is significant magnetoe- 
lastic coupling. It is therefore striking that small disor- 
der levels (x > 0.03) give rise to spin glass order at low 
temperature in place of the Neel phase. Disorder in the 
strength of exchange interactions, induced by distortions 
generated around Cd siteSfi^ seems a likely origin for this 
behaviour. 

Against this background, our aim in this paper is to 
study spin glass ordering in geometrically frustrated an- 
tiferromagnets with weak exchange randomness. Earlier 
studies of model frustrated systems in which low levels of 
disorder induce a spin glass phase are reported in Refs.fl^ 
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and[T^, and an earlier investigation by others of the prob- 
lem we consider is described in Ref. [la- A short account 
of some of our work has been given in a previous publica- 
tionii by two of the authors. Here we present extended 
results, including a mapping to an analogue of the con- 
ventional spin glass theory. We show, both analytically 
and from simulations, that weak exchange randomness 
indeed generates spin glass order, with a transition tem- 
perature proportional to the amplitude of exchange ran- 
domness, albeit with a different proportionality constant 
than in a conventional system without the geometrical 
frustration. As a result of the dominant, average antifer- 
romagnetic exchange, thermal fluctuations near the spin 
glass transition temperature are highly constrained: they 
lie within the ground state manifold of the equivalent sys- 
tem without quenched disorder. The ground states of the 
disorder-free system are macroscopically degenerate and 
have power-law correlations. Restriction of fluctuations 
to this set of states has two consequences. First, random- 
ness in nearest neighbour exchange generates effective in- 
teractions that are long ranged, though not sufficiently 
so to change the universality class for critical behaviour. 
Second, the transition temperature arising from a given 
amplitude of disorder is higher in a geometrically frus- 
trated system than it would be if the mean exchange 
interaction were zero. 

An outline of the paper is as follows: in Section [II] 
we introduce and discuss the replica treatment of a 
frustrated antiferromagnet with weak, homogeneous ex- 
change randomness; in Section IIIII we consider dilute 
disorder, for which we map the geometrically frustrated 
magnet to a set of pseudospins with random dipolar inter- 
actions; in Section llVI we present results from numerical 
simulations of the Heisenberg antiferromagnet on the py- 
rochlore lattice with exchange randomness; in Section IVl 
we use finite size scaling to analyse these results; and in 
Section [VII we make concluding remarks. 

II. REPLICA THEORY 

A. Setting the stage 

Our starting point is a pure system consisting of clas- 
sical m-component spins on the pyrochlore lattice with 
antifcrromagnctic interactions, described by the Hamil- 
tonian "Ho = J^i'^j- Here the sum T^^jj^ runs over 

pairs of nearest neighbours, and spins are of unit mag- 
nitude. This model is geometrically frustrated and has 
a macroscopically degenerate ground state. As a conse- 
quence, ordering is suppressed and the system remains 
paramagnetic down to zero tempcraturoi^ (except for 
the case m = 2, in which thermal fluctuations induce 
collinear order at low temperatures and which we ex- 
clude in the following). Another feature of the model is 
the emergence of power-law spin correlations at tempera- 
tures T Jfi^ and one of our concerns is to understand 
how these correlations influence spin glass phenomena in 



frustrated magnets. 

We introduce exchange disorder by replacing J in Hq 
with Jij = J + 5Jij , where dJij is random and has zero 
mean. We take the distribution of 6 Jij to be Gaussian 
with variance in our analytic work, and uniform on 
[—A, A] in our numerical simulations, setting A = A/VS 
so that the variance is the same in both cases. The regime 
of interest is A ^ J; the opposite situation, A ^ J, 
describes a system where disorder dominates over the 
antifcrromagnctic coupling and a conventional spin glass 
phase is expected at low temperatures. 

Wc hence consider the Hamiltonian 

H = J ^ S, • S, + ^ 6J,,S, ■ S, . (1) 

This model is of direct physical interest in both the 
Heisenberg and Ising cases, the latter arising as a de- 
scription of spin ice materials in which strong single ion 
anisotropy constrains spin orientations. 

B. Spin ice 

We first consider Ising spins, for which S; = ±1; the 
vector case is treated in Sec. Ill CI We use the replica 
tricki^ to carry out the disorder average. This produces 
an effective replica Hamiltonian 

1 A2 
RU — \ " qaca _ \ " qaobcacf) (r)\ 

{ij),a (ij),{ab) 

where the replica labels take values a,b = 1, ...n, the 
summation X](o&) over a <h, and the limit 7i — > is to 
be taken at the end of the computation. The last term 
can be rewritten as 

A^ 

ij,(ab) 

where Kij is the adjacency matrix on the pyrochlore lat- 
tice. The factor 1/2 accounts for double counting in sum- 
mation over sites. 

Next we introduce local overlaps defined on sites of 
the lattice via a Hubbard-Stratonovitch transformation^^ 
which decouples the four-spin term, yielding2i 

/3^Q = E U E QabK-j'Qi, - ^ E QafcS.-S? I . 

(ab) y ij * J 

(4) 

The overall partition function can now be written as 
an integral on of cxp[— /SHcff], where the effective 
Hamiltonian is given by 

/3^eff -\T. QabKr/Qib - [Q] ' (5) 

(ab)Aj 
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with 

M^^[Q] =log^eTE,„.),QL.S?S^^^ . (6) 

Here (• • • )o denotes an average with respect to the clean 
system Hamihonian, Hq. To advance further we need 
to evaluate VF^[(5]. An exact calculation is not possi- 
ble and so we rely on a high-temperature expansion of 
That is, we expand the exponential in powers of 

A/T. Since Hq does not couple different replicas, av- 
erages (••■)o with distinct replica labels factorise into 



It is convenient to split (3 Hon into S2, the part 
quadratic in Q, and Sint, an interaction part. Fourier 
transformation block diagonalises 82- Since the py- 
rochlore lattice has fours sites per unit cell, Fourier com- 
ponents carry both a wavevector q and a sublattice label 
K. We use a hat to denote the Fourier transform of a 



S2[Q] 

and SintiQ] 
where V^„,tu(q, p) 



A phase transition is signalled by the vanishing at 
the mean field critical temperature of one of the 

eigenvalues of the kernel of 82- Physically, one expects 
the critical mode to have wavevector zero. Because 
of translational invariance, a spatially uniform vector 
Q^jj = ipab • is an eigenvector of the kernel at all tem- 
peratures. This mode is also separately an eigenvector of 
K(0), with an eigenvalue that we denote by A4(0), and of 

6*^^(0). The temperature is therefore determined 

by the equation 

In the regime of most interest, A ^ J, the temper- 
ature scale J below which dipolar spin correlations de- 
velop is much larger than the spin glass transition tem- 



products of averages with the same replica labels. Aver- 
ages over odd numbers of spins vanish due to the absence 
of ordering in the pure system. As a consequence there 
is no term linear in Q. At order (A/T)^ and (A/T)^ 
only two-spin correlators Gij = (SiSj)o appear in the ex- 
pansion. At third order this is due to the fact that bare 
correlators with an odd number of spin vanish, and those 
with an even number of spins greater than two are ex- 
cluded due to constraints on replica indices that prohibit 
having four or all replica labels equal. Truncating the 
expansion at third order, we obtain 



(7) 

I 

function, and write the spin-spin correlator of the pure 
system as Gi^jj{q) — (S„(— q)Sr,(q))o- We also introduce 
the shorthand G^^ (q) = /p G'„,,(q— p)Gk,,(p). With this 
notation we have 



(8) 

(9) 
(10) 

I 

perature. For temperatures T ^ yMF gpj^ fluctuations 
are effectively confined to the groundstate manifold of 
the disorder-free system, and disorder acts as a pertur- 
bation. Under these conditions we can use the results of 
Rcfs.lT8land[23lto write the two-spin correlation function 
approximately in terms of the normalised eigenvectors 
{^^5=1,2,3,4} of the adjacency matrix K. This matrix has 
two flat bands, with eigenvalues X^^i^2 = —2 that are in- 
dependent of wavevector, and two dispersive bands, with 
eigenvalues A3(q) < A4(q). Wc hav o^^'^^ 

G.,(q) = 2 ^ (C/5)„(q)(C/5)^(q) . (12) 

5=1,2 

This approximation and the orthonomality of the 
eigenvectors {C/J yields Y.^^ G2^(0) 8. Now A4(0) = 6 
(see Ref.lH), and so we obtain from Eq. ([TT|) the result 



WaIQ] = ^ E QabGlQi, + |J E V^MbQicQL , where V,,, = G,,G,,G, 

{ab),ij {abc),ijk 



^ E / ^^.(-q) 



{ab),K,r] *^ 

A3 

r3 



{abc) ,KrjLO 

GK,,(q)GKcu(p)G,,^(-q - p) • 
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rpMF ^ A\/T2. A conventional mean-field result for T^^ 
is also eontained in Eq. dTTl) when one sets J = 0. In 
this ease G^riiq) = S^,j, whieh implies ^"^,(0) = 4 
and T^^/^o = AV6. Thus, within the fr amework of our 
calculation, there is an increase in the value of by a 
factor of due to the correlations arising from uniform 
frustrationi^ The increase in the spin glass phase tran- 
sition temperature compared with that of a conventional 
system with J = has a simple physical explanation: the 
geometrical frustration severely reduces the phase space 
available for spin-glass-destabilising thermal fluctuations 
in the low temperature manifold. 

To describe the transition we retain only the branch of 
soft modes. Wc denote the corresponding eigenvector of 
the kernel of 5*2 by w" (q), and the associated eigenvalue 
by E{q). Writing (5af,(q) ~ '^ab{<i)v''{q), the effective 
theory for the spin glass transition is 

= / ^(q)l^abq)p- (13) 

(ab) 

[^^f] E / ^(q'P)</'a&(q)v'ac(p)<y3&c(-q-p), 



Crucially, the power law decay G(r) 



at large r is 



where 



y (q, p)^J2 ^-n- P)^" (q)^' (P)^" (-P - q) (14) 



In conventional spin-glasses E{c\) ~ E{Q)-\-Ac^ as q ^ 0, 
with A a positive constant. If this behaviour persists in 
our case, the conventional effective theory is retrieved. 
Direct diagonalisation of 5*2 for q 7^ is complicated and 
we rely instead on non-degenerate perturbation theory 
to study the small q behaviour of E{q). To first order in 
(A/r)^ the eigenvalue E{q) is given by 



1 



A4(q) 



^\ Y.Um I G.,(q-p)G,,(q)[/4^q) 

/ Kl) "'P 



We consider in turn the contributions to E{q)—E{Q) from 
each term in this expression. First, using an explicit form 
for A4(q)^ it is easy to check that A4 ^(q) — A4"^(0) ~ 
Aoq^ as q — >■ with > 0. Next, direct evaluation of 

(q) — [74(0)] shows that it vanishes faster than 
as q — 0. Analysis of the contribution from the second 
term therefore reduces to the evaluation of 



'[jj 'io.qp^j\^c.ypG^^{p)\G^.r,{v). (15) 

where Vq, = d/dpa- The integral on p is most easily 
understood in real space, where it takes the form 



fast enough that the sum in Eq. (jl6p converges, giving a 
finite result for a — j3 and zero by symmetry for a ^ (5. 
Moreover, the signs of Eqs. (|15p and (|16p combine to give 
an overall positive sign to the coefficient of q^. Taken 
together these results ensure that E{c\) ^ E{Q) + Ac^ 
with positive A. 

Within the context of the effective theory it is also 
justified to replace V{(\, p) by its limit for q, p — > 0. 
Taking this limit in Eq. (jl4p and using the asymptotic 
form for G we find 

F(q,p) = (P,k)' + (q,k)2 + (p.q)^ - (p, q)(p, k)(q, k) 

where p, q, and k denote unit vectors in the directions 
of p, q and k = p -I- q, respectively. 

Summarising these results, we have obtained an essen- 
tially conventional replica theory at mean field level. A 
critical theory, almost identical to the one derived for 
Edwards-Anderson model2^, also follows. The only dif- 
ference to the standard version is in the form of the inter- 
action term, which has a non-trivial wave-vector depen- 
dence in our case, originating from dipolar correlations 
present in pure system. The effective critical theory is 



S[ip]=Y, j {Ci'+T)\^al>?{ci)+ (17) 

/ ^(q, P) E ¥'ab(q)'^ac(p)¥'6c(-P - q) , 
J ap ^ , ^ 



[abc) 



where r ~ 1 - {T^^^ /Tf and g - {l^^/Tf^ f 



(16) 



C. Vector spins 

In this section we sketch the generalisation of the re- 
sults derived above for the generic case of m-component 
spins. The derivation follows closely that for the Ising 
case with minor modifications due to additional spin 
component labels. Local overlaps («) now carry an 
additional pair of spin component labels a, /? and sum- 
mation over ab is unrestricted. Also, at variance with 
the Ising case diagonal terms (a ~ h) now give non- 
constant contribution to the effective Hamiltonian and 
cannot be dropped. As a consequence one has to de- 
fine diagonal overlaps Q^f (i), which contribute to the 
expansion of [Q] . This in turn leads to appearance of 
four and higher spin correlators already at order three in 
the expansion of W^. However it affects only the form 
of the cubic term, and not that of the quadratic term 
which determines the value of critical temperature. The 
mean-field critical temperature expression generalises to 
Tj^F = AyT2/m. The critical theory is very similar to 



the spin ice case, being 



aa ab.aP 1 




To conclude, critical replica theories for disordered 
frustrated magnets at all values of m (except, as indi- 
cated above, m = 2) coincide with the conventional crit- 
ical theory for spin-giassesi^ The underlying geometri- 
cal frustration reveals itself in a remaining wave-vector 
dependence of the interaction vertex and an increased 
mean-field value of the critical temperature as compared 
to the conventional case. 



III. DILUTE IMPURITIES 

It is also interesting to consider a model of dilute dis- 
order, in which a low density of isolated tetrahedra have 
exchange interactions that are different for different pairs 
of spins within the tetrahedron. We show in this sec- 
tion that each such tetrahedron has a pseudospin degree 
of freedom. We find that entropic interactions between 
these pseudospins, mediated by spins in the remaining 
tetrahedra, have a dipolar form. In these way we arrive 
at a similar conclusion to the one reached in Section HIl 
but for a different version of the problem and by a differ- 
ent route. 

As a first step, consider a single tetrahedron taken from 
this lattice, with spins Si . . . S4 at the vertices. With 
all exchange interactions equal, its ground states are the 
configurations for which Si = 0. The spin stiffness is 
zero in this toy problem in the sense that, within the set 
of ground states, the orientations of a pair of spins can 
be chosen arbitrarily. The consequences of fiuctuations 
in 5Jij with amplitude A have been set out in Ref . [igI and 
[25I : generically, a unique ground state is selected (up to 
global spin rotations) in which all four spins are collinear 
and the total spin of the tetrahedron is zero. In such a 
configuration the four spins can be grouped into two fer- 
romagnctically aligned pairs, and energy is minimised by 
picking these pairs appropriately, as illustrated in Fig. [T] 
Non-zero 5Jij hence induce a ground state stiffness, since 
changes in the relative orientation of a pair of spins cost 
an energy C'(A). 

Extending this discussion, consider a pyrochlorc lat- 
tice in which a randomly selected, dilute subset of special 
tetrahedra have interactions of unequal strength, while in 
the remainder all 5Jij are zero. Provided dilution is suf- 
ficiently high, ground states are macroscopically degen- 
erate, and in ground states each special tetrahedron has 
collinear spins at its vertices. The orientations of these 
quartets of collinear spins at different special tetrahedra 
are independent, and constitute some of the ground state 
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Figure 1: A ground state configuration of spins at vertices of 
an isolated tetrahedron in which there are antiferromagnetic 
exchange interactions of strength J between pairs joined by 
solids lines, and of strength J— 5 J (with < (5J <C J) between 
pairs joined by dashed lines. The vector yia is perpendicular 
to both of the links between ferromagnetically aligned pairs 
of spins 



degrees of freedom. We label the special tetrahedra by 
A and specify these orientations with unit vectors a a in 
spin space. In addition, to characterise the realisation of 
quenched disorder we introduce unit vectors in real 
space, defined to be perpendicular to both of the links 
on the lattice that join ferromagnetically aligned pairs of 
spins in tetrahedron A, as shown in Fig. [T] (Since in a 
serves to define an axis. ±n^ arc equivalent.) These vec- 
tors are each aligned along one of the cubic axes of the 
pyrochlore lattice and arc quenched random variables. 
At temperatures T ^ A only the xva-, and not the mag- 
nitudes of the i5Jij, are important to characterise the dis- 
order. 

Integration over all other degrees of freedom induces 
an effective, entropic interaction between the cr^'s. The 
form of this effective interaction can be deduced by us- 
ing the parameterisation of the ground states of the 
disorder-free model in terms of a gauge field, introduced 
in Refs. Hi. 

The essence of this parameterisation can be sum- 
marised as follows, treating in the first instance the case 
of Ising spins Si. A three-component vector field B(r) 
is chosen to represent spin configurations, in such a way 
that the condition for a configuration to be a ground 
state is V • B(r) = 0. At the microscopic level, this is 
achieved in two stepsJ^ First, one notes that the cen- 
tres of tetrahedra on the pyrochlorc lattice themselves 
form a diamond lattice, which is bipartite. It is there- 
fore possible to define a unit vector e(vi) in real space 
at each site r.^ of the pyrochlore lattice, with the ori- 
entation convention that it is directed from one chosen 
diamond sublattice towards the other. Second, one de- 
fines B(ri) = e{ri)Si, which has zero lattice divergence in 
ground states. After coarse-graining, B(r) is treated as 
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continuous, divergence-free field. Extending these ideas 
to treat m-component spins, one introduces m fields 
B"(r). with a = 1 . . .to, related at the lattice level to 
spin components 5f by B"(ri) = e(yi)Sf. The coarse- 
graining procedure gives rise to an entropic weight that 
favours configurations with small field strengths. Writing 
this weight as e~'^°, 5*0 is postulatedi^ to have the form 

= d'rE|B"(r)P (19) 

where k characterises an entropic stiffness, which is dis- 
tinct from the energetic stiffness that arises when 5 J^- is 
non-zero for all nearest neighbour pairs. 

In this language, the condition that the spins of a spe- 
cial tetrahedron located at a random position are 
coUinear with orientation a a translates into the condi- 
tion on the fields B"(r) that B"(rA) = a\wA- We im- 
pose these constraints by introducing three-component 
fields and using on the special tetrahedra 

^ (B«(rA) - alx^A) = ^ y d3</,^e"^S (B°(r)-.2n^) 

Integrating out the fields B"(r) and the variables cr^, we 
arrive at a weight e"'^"'"' for 0^ of the form 

= ^ E • M(r^ - rs) • 0S + E ^(^i) • (20) 

olAB a 

Here the 3x3 interaction matrix M(r) is dipolar, with 
elements [M(r)]fc; = (r^bu - 3rA;n)/47rr^. The one- 
body term V{if>A) is a function of if \ = J2ai^A ' ^a)^ 
and has the expansion V{(p'^) = C2f'^ + c^ip'^ ■ ■ ■ , where 
C2 = (2to)"^ and C4 = {Am'^[m + 2])^^. This effective 
model has site disorder, since the tetrahedra labelled by 
A and B are selected at random. By this means, we 
have arrived at a model of randomly located interacting 
dipoles as a description in the T — !• limit. We expect 
the model to have a classical, zero-temperature phase 
transition between a paramagnetic phase at low density 
and a spin glass phase at high density. 

At finite temperature thermal excitations generate a 
finite correlation length ^, which sets a maximum range 
for the interaction M(r). This correlation length diverges 
in the low temperature limit, exponentially in J/T for 
the Ising spins, and as the power law ^ ~ {J/TY^'^ for 
Heisenberg spins. On increasing temperature from zero 
in the spin glass phase, a transition to a paramagnet is 
expected, with a transition temperature set by the lower 
of two scales: one of these is the temperature at which ^ 
becomes comparable to the spacing between defect tetra- 
hedra; the other is the disorder strength A. Despite the 
dipolar form of interactions at distances shorter than ^ 
this transition is expected to be in the same universality 
class as with short-range exchangei2& 

In principle an analogous finite-temperature mapping 
could be made in terms of the model of Section II, by 
replacing the hard constraints ^"{rA) = cr^nA with 



soft weightings exp {-Ea{{J^})/T). These will be de- 
termined by the energetic costs Ea{{J^}) of local spin- 
configurations drawn from the ground state manifold of 
the clean system, where {Jij} denotes the set of {Jij} 
within tetrahedron A. This mapping would yield an ana- 
logue of Eq. ((20)) with a summation over all A but with 
extra quench-random local weightings, again leading 
to spin glass behaviour as expected for a site-disordered 
Ginzburg-Landau spin glass. 



IV. NUMERICAL SIMULATIONS 

In order to investigate our ideas further we turn to 
Monte Carlo simulations of the classical Heisenberg an- 
tifcrromagnct with nearest ncigbour exchange on the py- 
rochlorc lattice. Our focus is on the effects of weak ran- 
domness in the strength of exchange interactions. This 
section provides a complete description of work presented 
briefly in an earlier publicationJi 



A. Model and Method 

We take the exchange interaction between spins at 
neighbouring sites i and j to have strength J + 6 Jij with 
6 Jij an independent random variable for each bond, uni- 
formly distributed in the range [—A, A]. Our interest is 
in the limit A ^ J and our most extensive results are 
for A = O.IJ. 

We simulate a sample in the shape of a rhomboid that 
has edges parallel to the primitive basis vectors of the lat- 
tice, with periodic boundary conditions between opposite 
faces. System size is specified by the linear dimension L 
of the sample. The number of primitive unit cells in such 
a sample is = L'^ and the total of number of spins is 
A^s ~ 4L^. We present data for sample sizes in the range 
from L = 2 to L = 7. At each sample size it is necessary 
to average over different disorder realisations. We used 
10'^ realisations for L = 2 but found that 200 realisations 
are sufficient for L = 7. 

We employ parallel tcmpcring^l to ensure equilibra- 
tion of large systems at low temperature. In this ap- 
proach, one simulates A"array copies of the system simu- 
lancously. Each copy is at a different temperature, taken 
in a range from Tmin to Tmax with geometric spacing. 
This range is required to be wide, since Tmin must be 
below the spin glass transition temperature while Tmax 
must be high enough that relaxation at that temperature 
is fast. At the same time, adjacent temperatures should 
be sufficiently close that there is a high probability for 
configurations to be exchanged between them under the 
moves of the parallel tempering algorithm. This requires 
a sufficiently large value value of A'array • We take^i 

^array = K^^' ln(T,„ax/T„,„) . (21) 
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Parallel Tempering Steps 



Figure 2: Evolution of the spin glass correlation function 
for two different starting configurations. Simulations start- 
ing from a random initial condition (for r = 1 (A), r — rmax 
(o)) are compared with simulations starting in a Neel ordered 
state (for spatial separations r = 1 (o), r = rmax (d))- Error 
bars are omitted for clarity: in the worst case (r = )"max) they 
are ~ 10%. 



Figure 3: Cv against L for two temperatures: that of the 
maximum (o); and T/A = 0.1 (A), the lowest temperature 
simulated for larger system sizes. 



over them as in Ising systems. 



B. Testing Equilibration 



We estimate equilibration times by studying the evo- 
lution of observables starting from different initial states. 
Two simple choices of initial state are an infinite temper- 
ature configuration with random spin orientations, and a 
Neel ordered configuration with collinear spins, which is 
a ground state of the model without exchange random- 
ness. For each of these initial states we show in Fig. [2] 
the evolution with Monte Carlo time of the spin glass 
correlation function C(r) [defined in Eq. ([H])] for r = 1 
and r = rmax, the maximum separation in a sample of 
size L = 7, taking A/J = 0.1 and T/A = 0.1. This is 
the largest lattice size and lowest temperature, and hence 
the most difficult case, for which we present detailed re- 
sults in our study of the spin glass transition. As seen 
from Fig. [2l the equilibration time in this case is ~ 10^ 
parallel tempering steps, although memory of the differ- 
ence between initial configurations is lost after a shorter 
time. When using parallel tempering it is also important 
that each copy of the system should visit every tempera- 
ture simulated with equal probability. We have checked 
that the simulation time for a copy to loose memory of 
its initial temperature is shorter than the equilibration 
time for the spin glass correlation function. On the basis 
of these tests, for L = 7 we collect data after an equili- 
bration time of 10^ parallel tempering steps. For smaller 
system sizes equilibration is more rapid. 

Our equilibration times are similar to those for conven- 
tional Heisenberg spin glasses on the cubic lattice^. It 
is worth noting that these equilibration times are smaller 
than those typically needed for Ising spin glasses. Lee 
and Young suggest that this is because energy barriers 
are smaller in the Heisenberg modeli^ The extra degrees 
of freedom in the Heisenberg model mean that the simu- 
lation can find paths around energy barriers, rather than 



C. Specific Heat 

The low temperature heat capacity Cy of the classical 
Heisenberg antiferromagnet on the pyrochlore lattice is 
interesting as a diagnostic for macroscopic ground state 
degeneracyii^ Without degeneracy, equipartition and the 
fact that each spin has two degrees of freedom would give 
a classical low temperature heat capacity of A;b per spin. 
The smaller measured value of (3/4) fee per spin demon- 
strates that one quarter of the degrees of freedom in the 
model make no contribution to Cy because they can fluc- 
tuate without energy cost. We expect exchange random- 
ness to eliminate this macroscopic number of zero modes, 
leaving only the three zero modes associated with global 
spin rotations. The limiting low temperature value of 
the heat capacity per spin should then be (1 — 3/2A^s)A;B■ 
At higher temperatures the heat capacity is expected to 
have a broad maximum in the vicinity of the spin glass 
transition, and to remain finite and smooth even in the 
thermodynamic limit. 

In our simulations we determine the heat capacity 
from the variance of energy fluctuations. As reported 
previouslyjii the heat capacity has a broad maximum 
around T/A ^ 0.45 for A/J = 0.1, while at low tem- 
peratures it tends to unity for large system sizes. We 
show in Fig. [3] the dependence on system size of at 
the temperature of the maximum and at the low temper- 
ature T/A = 0.1. As expected, both values approach a 
constant with increasing system size, which in the sec- 
ond case is close to fee. Further calculations down to 
T/A = 10~^ for L = A (not shown) confirm the expected 
value (1 — 3/2Ns)kB more precisely. 
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Figure 4: C(rmaa;) versus temperature for L = 7(A), L = 
6(n), and L = 5(o). The inset shows the scaling of C{rmax) 
with system system size for T/A = 0T(o), T/A = 0.2(n), 
and T/A = 0.4(o). 



D. Spin Glass Correlation Function 

To search for spin freezing we study the spin glass cor- 
relation function C(r). This is defined in terms of be- 
haviour in two copies of the system with identical dis- 
order. Denoting thermal averages in copies 1 and 2 by 
(. . .)i and (. . .)2, and the disorder average by [. . .]av, we 
have 

C(r) = [(S(0)-S(r))i(S(0).S(r))2]„„. (22) 

Spin freezing is indicated by a non-zero limiting value 
for C(r) at large r. We show in Fig. |4] the tempera- 
ture dependence of C(r) for the maximum spin separa- 
tion (r ~ ^niax) hi the three largest system sizes studied 
{L = 5,6, 7). There is a clear transition within the tem- 
perature range 0.2 < T/A < 0.4. This behaviour is in 
marked contrast to that of the pure system, where C(r) 
falls with r as in the low temperature limit, and ex- 
ponentially at finite temperaturci^ 

The inset to Fig. S] suggests that C(rmax) tends to a fi- 
nite constant below Tj for large L, further supporting our 
conclusion that there is a finite temperature transition in 
the infinite system when weak bond disorder is present. 
The behaviour of the correlation function puts simple 
bounds on the transition temperature, 0.2 < T/ /A < 0.4. 

The behaviour of the spin glass correlation function 
thus provides convincing evidence for a spin glass tran- 
sition in the model. However, the value of the critical 
temperature has a large uncertainty. As we discuss in 
the next section, finite size scaling can be used to sharpen 
the estimate for the critical temperature and to find ap- 
proximate values for critical exponents. 

V. FINITE-SIZE SCALING 

There is extensive past work on finite size scaling anal- 
ysis of the transition for conventional models of spin 



Figure 5: The spin glass susceptibility, XSG vs T/A for system 
sizes from L = 4 (closed triangles) to T = 7 (open circles). 



glasses that have zero mean exchange interactioni^ From 
this it has emerged that study of the behaviour of the 
spin glass correlation length is a particularly effective ap- 
proach. We find that the situation is different in the case 
of geometrically frustrated systems. These have strong 
short-range correlations at low temperature, even in the 
absence of disorder, which complicate a scaling analysis 
using the limited range of system sizes available. Instead 
we employ scaling collapse of the spin glass susceptibility 
[Eq. ([23|] to obtain simultaneous estimates of Tf, i' and 
7. We emphasise that our aim here is not to find precise 
values for the critical exponents; rather, our principal 
objective is to confirm that there is a finite temperature 
second order phase transition and determine its temper- 
ature. 

The spin glass susceptibility, which is related to the 
non-linear susceptibility Xni-, is defined as^ 

XSG-E^W- (23) 

r 

In the paramagnetic phase XSG ^ C'(l) since the only 
significant contribution to C(r) is from small r. Its di- 
vergence at finite temperature signals a phase transition. 
Approaching the critical temperature from above, we ex- 
pect 

XSG - t-^ , (24) 

where t ^ [T — Tf)/Tf and 7 is the corresponding crit- 
ical exponent. Eq. (IMl) holds if the dimension d of the 
system is greater than d;, the lower critical dimension, 
so that T{ > 0. The evidence from the previous section 
strongly favours this scenario. Furthermore, from the re- 
sults of Section |TT] we expect the upper critical dimension 
for the transition to have its conventional value, which 
is believed to be sb^. The hyperscaling relations should 
therefore hold, and using di/ = 7 -I- 2/3 we can deduce the 
value of /3 from scaling of xso- 

If there is a thermodynamic phase transition then we 
anticipate close to Tf in a finite system the scaling be- 
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Figure 6: Scaling collapse of XSG for: (top panel) v = 1.2, 7 
1.2,Tc/A = 0.2; (middle panel) u = 1.0,7 = 1-45, Tc/A 
0.23; (bottom panel) v = 0.9, 7 = 1.6, T^/A = 0.32 



haviour 



x(r,i) = L^/V(i'/''t) 



(25) 



where v is the critical exponent of the correlation length. 
The scaling function obeys /(O) > and f{x) ~ x^'^ for 
X — >■ cxj. In Fig. Owe show as a function of T/A for 
system sizes from L = 4 to L = 7. 

The rapid increase in xsG below T/A ~ 0.4 provides 
clear evidence of a spin glass transition. The scaling anal- 
ysis of xsG is complicated by large finite-size effects. Due 
to the ground state constraint in frustrated antifcrromag- 
nets, our model has significant short-ranged correlations 
that are approximately independent of system size. For 
L ~ 2 and L ~ 3 the contribution to xsG from these local 
correlations is significantly greater than the contribution 
from C(r) at large distances and so these system sizes 
cannot be included in the analysis. 

A scaling collapse of xsG for the system sizes 4 < L < 7 
yields the best-fit parameters Tf/A = 0.23(9), v = 1.0(2) 
and 7 = 1.45(45)^. It is difficult to make detailed esti- 
mates of the errors in these values, but a simple approach 
is to explore the range of parameters that still gives rea- 
sonable data collapse. In Fig. [6] the scaling collapse for 
the best fit is compared with behaviour for two 'worst- 
case' fits: (i) Tf/A = 0.2, = 1.2 and 7 = 1.2; and (ii) 
Tf/A = 0.32, V = 0.9 and 7 = 1.6. These worst-case 
fits were generated by fixing the value of Tf and then 
adjusting the values of u and 7 to minimise scatter of 
the data. As the data collapse is visibly poorer for both 
the worst cases, we believe they set bounds on the value 
of Tf. In addition, since the exponent values that pro- 
duce good scaling collapse are correlated with the fitting 
value for Tf , they give bounds on v and 7. These are the 
uncertainties quoted above. 

Our results can be compared with those from simula- 
tions on conventional spin glass models, and with exper- 
iment. As a first step, it is interesting to examine the 
value of Tf. The most direct comparison would be be- 
tween Tf/A at large J and at J = 0, evaluated in both 
cases for the Hcisenbcrg model on the pyrochlorc lattice. 



Because we do not have data for this lattice at J = 0, we 
compare instead with established results for the Heisen- 
berg spin glass on the simple cubic lattice. Since the 
number of nearest neighbours is the same on both lat- 
tices, we expect that this comparison will be adequate to 
establish the trend in Tf/A with J. For the cubic lat- 
tice with Gaussian nearest neighbour exchange of zero 
mean and unit variance, Tf = 0.129,^ Assuming that it 
is appropriate to compare our rectangular exchange dis- 
tribution with a Gaussian by simply equating variances, 
we conclude that large J increases the value of Tf/A by 
the factor V3 x 0.23/0.129 « 3.1. This substantial effect 
is physically reasonable: thermal fluctuations in a ge- 
ometrically frustrated system are much more restricted 
than in a conventional spin glass. The spin glass phase 
therefore extends to higher temperatures than in a sys- 
tem with average exchange J = 0. Turning to exponent 
values, our result for ly lies within the range (1.01 — 1.50) 
reported from simulations of the conventional Hcisenbcrg 
spin glass on the simple cubic latticei^ Gomparison with 
experimental results for the pyrochlore antiferromagnet 
Y2M02O7 is also possible. The experimental values are 
7 = 2.9(5) and f3 = 0.8(2)^. Our value of 7 (given above) 
is significantly smaller than the experimental one, while 
the result wc obtain for /3 using hypcrscaling, /3 = 0.8(3), 
is in close agreement. Experience from very large scale 
simulations of conventional three-dimensional Heisenberg 
spin glasses unfortunately suggests that a much bigger 
computational effort than our own would be necessary 
to characterise critical behaviour reliably. We neverthe- 
less have confidence in our central result from simula- 
tion: that weak exchange randomness in a geometrically 
frustrated magnet induces a spin glass transition, with a 
higher transition temperature than would be the case of 
the mean exchange were not strongly antiferromagnetic. 



VI. CONCLUSIONS 

Spin freezing has been observed experimentally at low 
temperatures in many geometrically frustrated magnetic 
materials, and its origin has long been unclear. In this 
work we have carried out a detailed analytical and nu- 
merical study of models with weak disorder, showing that 
this is a mechanism that produces freezing and a low tem- 
perature spin glass phase. 

Analytically, we have studied geometrically frustrated 
antiferromagnets perturbed by weak exchange random- 
ness. This form of disorder may be generated by random 
strains in the sample, arising from non-magnetic chemi- 
cal disorder. Wc have shown that it leads to a spin glass 
phase at low temperatures. The main result, expressed 
by Eqns. ((TT)) and (jl8p . is that a model of the Edwards- 
Anderson (EA) type, but with dominant, mean geomet- 
rically frustrated antiferromagnetic exchange, falls into 
the same universality class as the usual EA modeliiS 
We find no essential deviations from conventional spin 
glass results for the character of the transition to the low 
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temperature phase. However, the (mean-field) transition 
temperature is increased by the long-range correlations of 
the pure antiferromagnet, by a factor of \f2 as compared 
to a system with no mean antiferromagnetic interaction. 

In a complementary approach, we have shown quali- 
tatively how a similar conclusion arises in a model with 
dilute disorder. For this case we map the initial system, 
consisting of spins on a regular, frustrated lattice with 
a low concentration of disordered interactions, to an ef- 
fective system, made up of pseudospins at random sites 
with entropic interactions that arc dipolar in character. 
Such interactions, although long range, are expected to 
yield a conventional spin glass transition;^ as we find for 
homogeneous disorder. 

In Sec llVl we have checked these ideas using numerical 
simulations. For conventional Ising Edwards- Anderson 
systems it has long been accepted that there is a finite 
temperature spin glass transition in three dimensions. 
There is also compelling evidence for spin glass ordering 
in real experimental three-dimensional spin glasses. The 
existence of a spin glass transition in three-dimensional 
Heisenberg Edwards- Anderson spin glasses has been con- 
troversial, but the most recent exhaustive studies^ in- 
dicate that such a transition does occur. The data we 
present supports the corresponding conclusion that spin 
freezing transitions also occur in disordered geometrically 
frustrated antiferromagnets. A disputed issue for three- 
dimensional Heisenberg Edwards- Anderson spin glasses 
is whether the spin-glass transition is distinct from a 
chiral ordering transitioni^ We have not attempted to 
address this for geometrically frustrated systems but we 
would not expect any qualitative differences from con- 



ventional spin glasses without geometrical frustration. 

While our simulations have been of a three-dimensional 
system and the replica treatment we have presented 
has been at a mean-field level, some geometrically frus- 
trated antiferromagnets exhibiting spin freezing, includ- 
ing SrCr8Ga40i9, are in fact quasi- two dimensional. 
Since our central conclusion is that universal features 
of spin glass ordering in geometrically frustrated mag- 
nets with weak quenched disorder should be the same 
as those in conventional Edwards- Anderson systems, and 
since the two-dimensional Heisenberg Edwards- Anderson 
model is believed not to have spin-glass ordering^ it is 
necessary to appeal to weak interlayer coupling or spin 
anisotropy to account for spin freezing in quasi-two di- 
mensional geometrically frustrated Heisenberg systems; 
note that the relevant scale for the interlayer coupling to 
affect spin glass ordering is set by A rather than J. 

Our numerical simulations support the theoretical ex- 
pectation that the critical temperature is proportional to 
A for J > A, with T/ ~ 0.23A ~ 0.40A. The observed 
values of T/ therefore imply fluctuations in exchange in- 
teractions strength with a variance that is a few percent 
of the mean in SrCr8Ga40ig and some ten times larger 
in Y2M02O7. A direct experimental search for such ex- 
change fluctuations would be of great interest. 
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